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We study spectral singularities and their application in determining the threshold gain coefficient 
g ( E / M ) f or oblique transverse electric/magnetic (TE/TM) modes of an infinite planar slab of ho¬ 
mogenous optically active material. We show that g ^ is a monotonically decreasing function of 
the incidence angle 8 (measured with respect to the normal direction to the slab), while has 
a single maximum, 8 C , where it takes an extremely large value. We identify 9 C with the Brewster’s 
angle and show that and g l ' M ' 1 coincide for 6 = 0 (normal incidence), tend to zero as 6 —► 90°, 

and satisfy g ^ < g for 0 < 8 < 90°. We therefore conclude that lasing and coherent perfect 
absorption are always more difficult to achieve for the oblique TM waves and that they are virtually 
impossible for the TM waves with 8 ss 8 C . We also give a detailed description of the behavior of 
the energy density and the Poynting vector for spectrally singular oblique TE and TM waves. This 
provides an explicit demonstration of the parity-invariance of these waves and shows that the energy 
density of a spectrally singular TM wave with 8 > 8 C is smaller inside the gain region than outside 
it. The converse is true for the TM waves with 8 < 9 C and all TE waves. 

Pacs numbers: 03.65.Nk, 42.25.Bs, 42.60.Da, 24.30.Gd 


I. INTRODUCTION 

The history of modern physics includes happy instances where a purely mathematical concept turns out to describe 
an interesting physical phenomenon. A very recent example of such an instance is the remarkable observation that 
the spectral singularities of non-Hermitian differential operators [l|, which constituted a subject of study in pure 
mathematics for over half a century turn out to find concrete physical realizations in terms of lasing at the 
threshold gain Q and coherent perfect absorption (CPA) of electromagnetic waves Q. 

Spectral singularities are certain points of the continuous spectrum of non-Hermitian operators. They draw the 
attention of mathematicians because they are responsible for a number of mathematical peculiarities that can never 
arise for Hermitian operators. In the context of the scattering theory, spectral singularities correspond to scattering 
states (with real and positive energy) that behave exactly like resonance states; they are zero-width resonances 
[3|. This observation has provided ample motivation for the study of their physical aspects fsl lTs^ and nonlinear 
generalizations fl9l - [23j . For a recent review, see J24| . 

An interesting outcome of the study of spectral singularities in optics is a mathematical derivation of the lasing 
threshold condition for a variety of optical setups. Among the notable applications of this approach is the determi¬ 
nation of the threshold condition for the radial modes of single-layer and bilayer material with spherical geometry 
EMU, the TE modes of 'PT-symmetric and non-PT-symmetric bilayer slabs EHHHj and the whispering gallery 
modes of cylindrical and spherical gain media 0,0. 

The fact that the mathematical condition for the emergence of a spectral singularity implies the laser threshold 
condition was initially shown in an explicit manner in |6[| for the TE modes of an infinite planar slab laser where the 
wavevector k is taken along the normal direction e z to the boundary of the slab. A natural extension of the approach 
of @ is to consider the more general case where the direction of k deviations from e 2 , i.e., investigate the behavior of 
oblique waves (See Fig. [T]) Ref. EB treats this problem for the TE modes. The purpose of the present paper is to 
offer a more systematic and general solution of this problem that applies for both TE and TM modes. In particular, 
we show that much of the analysis of EH can be reduced to that of Q (i.e., to the case of normal incidence) by a 
simple change of variables, derive an explicit form of the laser threshold (and CPA) condition for oblique TE and 
TM modes, and examine the behavior of the Poynting vector and the energy density of the waves for the spectrally 
singular configurations of the system. This reveals a peculiar property of the spectrally singular TM waves having 
an emission angle larger than the Brewster’s angle; the magnitude of the Poynting vector and the energy density of 
these waves are smaller in the gain region than in the surrounding vacuum. 
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FIG. 1: (Color online) TE (on the left) and TM (on the right) modes of an infinite planar slab of thickness L. 



II. TE AND TM MODES OF AN INFINITE PLANAR SLAB 


Consider an infinite planar slab of thickness L that is aligned in the x-y plane of some Cartesian coordinate system 
{(x,y,z)}. Let n denote the complex refractive index of the content of the slab, and suppose that it is homogeneous 
and time-independent. Maxwell’s equations describing the interaction of the electromagnetic waves with this system 
have the form: 

V -V = 0, 
d t V - V x U = 0, 

where £ and T~L are the electric and magnetic fields, 


o' 

II 

f £Q 

r> 

(1) 

d t B + V x £ = 0, 

(2) 


£o and fi o are respectively the permeability and permittivity of the vacuum, and 


f n 2 for 0 < z < L, 

| 1 for z < 0 and z > L. 


( 3 ) 


For time-harmonic electromagnetic fields, £(r,t) = e lbJt E(r) and %(r, t) = e luJt H(r ), Maxwell’s equations reduce 
to 

[V 2 + k 2 j(z)] E(r) = 0, H(r) = -^-VxE(r), (4) 

[V 2 + k 2 i(z)\ H(f) = 0, E(r) = x H(f), (5) 

where r := ( x : y,z ), k := io/c represents the wavenumber, and c := 1/y/JIoeo and Zq := \Zyo/eo are respectively the 

speed of light and impedance in vacuum. 

The transverse electric (TE) and transverse magnetic (TM) waves correspond to the solutions of (|4]) and J5]) for 
which E[f) and H(r) are respectively parallel to the surface of the slab. Let us choose our y-axis to be along this 
direction, and demand that outside the slab E(f) (respectively H(f)) coincides with a plane wave with wavevector k 
in the x-z plane, i.e., 

k = k x e x + k z e z , k x := ksinO , k z := k cos#, (6) 

where e x ,e y , and e z , are respectively the unit vectors along the a>, y- and z-axes, and 8 € [—90°,90°] is the incidence 
angle (See Fig. [T]) Then, the electric field for the TE waves and the magnetic field for the TM waves are respectively 
given by 

E{f) = £’(z)e ik * x e y , H(r) = Jt°(z)e ik * x e y , 


( 7 ) 
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where $ and 3^ are solutions of the Schrodinger equation 

-i/)"(z) + v(z)ip(z) = k 2 ip(z), z£{ 0,L}, (8) 

for the potential v(z) := k 2 [ 1 + sin 2 9 — 3 ( 2 ;)]. Because v{z) is a piecewise constant potential, we can easily solve © 
to obtain 


1 

r a 0 e ikzZ + 6 0 e~ lkzZ 

for 

z < 0 , 


'ip(z) := i 

ai e ikz + bi e~ ikz 

for 

>4 

V 

V 

0 

(9) 

1 

{ a 2 e lkzZ + b 2 e~ lkzZ 

for 

z > L, 



where at and bi , with i = 0 , 1 , 2 , are complex coefficients, and 

k := k\/n 2 — sin 2 9 = fc 2 n, n := ^ ~ . (10) 


In particular, $(z) and J^(z) are given by the right-hand side of © with generally different choices for the constants 
di and bi. 

Substituting (© in the second equation in (© and ©, we can find the magnetic field for the TE waves and the 
electric field for the TM waves. The resulting expressions are acceptable provided that they satisfy the appropriate 
boundary conditions for the problem, namely that at z = 0 and z = L the tangential components of E and H must 
be continuous functions of z. Table U gives explicit expressions for the components of the electric and magnetic fields, 
and Table nil lists the corresponding boundary conditions. They involve the following quantities. 




ST[x,z) 


u 


r °o e ikzZ — 

< ai e ikz - 
{ a 2 e ikzZ - 


b Q e~ ikzZ for 
bi e~ ikz for 
b 2 e~ ikzZ for 


z < 0, 

0 < z < L, 
z > L. 


\/n 2 — sin 2 9e ikxX for z€ [0,L\, 
cos 9e ik * x for z ^ [0, L], 



V n 2 — sin 2 9 

for TE waves, 

n = 

cos 9 

ft 

\/n 2 — sin 2 9 

for TM waves. 



n 2 

n 2 cos 9 



( 11 ) 

( 12 ) 


(13) 


TE-Fields 

TM-Fields 

Ex — E z = Hy = 0 

Ey = £{z) e ikxX 

H x ^~^E-3T(x,z) 

TT sin 6e ikxX £(z) 

H “ Zo 

E y — H x = H z = 0 

_ Zosm9e ikxX J^(z) 

E ‘ ‘~ >w 

H y = Jf(z) e ikxX 


TABLE I: Components of the TE and TM fields in cartesian coordinates. Here <a(z) is given by the right-hand side of j9j, and 
J?( 2 ) and 3T{x,z) are respectively defined by (lllll and (1121) . 


III. TRANSFER MATRIX AND SPECTRAL SINGULARITIES 


The transfer matrix for the system we consider is the 2x2 matrix M satisfying 


a 2 

b 2 



( 14 ) 
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0 

II 

ao + bo = ai + bi, bo — ao = u(6i — a 1 ) 

II 

ikL , 7 —ikL ik z L , 7 —ik z L 

a ie +b ie = a 2 e + tee 

/ ikL 7 —ikL\ ik z L 7 —ik z L 

u(a ie — b\e ) = 


TABLE II: Boundary conditions for TE and TM waves. Here u is defined by Ga¬ 


in view of the fact that the (left/right) complex reflection and transmission amplitudes, R l ^ r and T l ^ r , are defined by 

Left-incident waves (62 = 0): 

Right-incident waves (ao = 0): R r := T l := 


R l := , 

Oo 

BT := 

02 


T l := —, 
ao 


Equation (TbO) relates the entries of M to R l ^ r and T l / r according to 


R l = — 


Mo 


R r = 


Mi 


rj~il _ 


det M 


M22 M22 M22 

With the help of the boundary conditions given in Table [Til we can easily show that 


rj-ir _ 


M 


22 


(15) 


M n = 

|cos(fcL) + ^(u + u 1 )sin(fcL)|e 

(16) 

M\2 = 

^(u — u -1 ) sin (kL)e~ lkzL , 

(17) 

M 2 i = 

— ^(u — u _1 ) sin (itL)e lkzL , 

(18) 

M 22 = 

|cos(/cL) — -(u + u _1 ) sin(fcL)| e lkzL . 

(19) 


In particular, we have 


det M = 1. 


( 20 ) 


For TE waves the transfer matrix M coincides with the one defined by the Schrodinger equation © and the require¬ 
ment that the solution be continuously differentiable at z = 0 and z = L. In this case (l20l) follows from well-known 
Wronskian identities of second order linear differential equations [25j, [2(|. These do not however apply for the TM 
waves, because TM waves are given by solutions of the Schrodinger equation © involving jump conditions at z = 0 
and z = L. We can identify the latter with the effect of certain point interactions at these points. It is known that the 
transfer matrix for a point interaction does not generally satisfy ®, 53- The point interactions giving rise to the 
TM jump conditions at z = 0 and z = L are separately examples of the anomalous point interactions which violate 
However, one can show (using the results of [13]) that the contribution of each of these point interactions to 
det M cancels the other’s, so that (l20l) holds. 

The spectral singularities correspond to the real and positive values of the wavenumber k for which M 22 = 0. This 
implies that the reflection and transmission amplitudes diverge, jlj. For the normally incident TE waves, where 


= 0, 


k , = k. 


k = fcn, 


u = n = n. 


( 21 ) 


Equations m © reproduce the known expression for the transfer matrix of a complex barrier potential, and 
M 22 = 0 takes the form [28| : 


-2ikLn 


n- 1 
n +1 


( 22 ) 


In view of (1191) . we can obtain the spectral singularities in the oblique TE and TM modes of our system by making 
the following substitution in 


k^-=kdl- 


sin 


n —> u = 


n for TE modes, 


— for TM modes. 


(23) 
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This gives 


g—2 ikL _ 


u- 1 



(24) 


Next, we substitute m and (fl3l) in (1241) to obtain the following more explicit relation for the spectral singularities 
in oblique TE and TM modes of our model. 


g—2 ikLyJ n 2 —sin 2 9 


\J n 2 — sin 2 9 — cos 9 
\/n 2 — sin 2 9 + n e cos 9 



£:= 


0 

2 


for TE modes, 
for TM modes. 


(25) 


The expression given in (1251) for the TE modes differs from the result reported in Ref. [Ts) . The root of this discrepancy 
is that the author of [18| takes different wave vectors inside and outside the slab and relates their direction using the 
usual expression for Snell’s law, namely nsind = n' sin (9'. Clearly, this cannot hold for complex refractive indices with 
different phases. We should also like to emphasize that a proper analysis of this problem should only rely on the use 
of Maxwell’s equations and the boundary conditions on the electromagnetic fields. In particular, there is no need for 
invoking the use of consequences of the latter such as Snell’s law and deal with the subtilises associated with their 
correct generalization and proper application. This is the basic strategy we pursue in the present article. 


IV. THRESHOLD GAIN 


By definition the gain coefficient is given by 


g := —2kn = — 


47 TK 

~ 


(26) 


where n is imaginary part of the refractive index n and A := 2n/k is the wavelength. The gain coefficient required for 
the emergence of a spectral singularity is known as the threshold gain coefficient @. For the oblique TE/TM modes 
we denote this quantity by and use (1251) to derive an expression for it. This is known as the lasing threshold 

condition in optics, |30j. 

Let us denote the real part of n by g , so that 


n = g + in. 


and take the logarithm of both side of (1251) to obtain 


i ( n' — n e cos 9 A 
2 Ln' \ n' + cos 9 J 


(27) 


(28) 


where n' := \/n 2 — sin 2 9. If we solve (1261) for k and substitute the result and (1271) in (1251) . we obtain a complex 
equation involving g , g, k, and 9. The left hand-side of this equation is real, so it should be equal to the real part of 
the right-hand side, and the imaginary part of the right-hand side should vanish. This gives a pair of real equations 
that can be manipulated to arrive at 


g (E/M) 


2 Im(n) j n' — n f cos(? 
L Im(n') n' + n e cos 9 
27rLRe(n') 

7 Tin — ip( E / M ) ’ 


(29) 

(30) 


where Re and Im respectively stand for the real and imaginary parts of their argument, X^ E / M ^ is the wavelength of 
the spectrally singular TE/TM wave generated for the gain coefficient g( E / M \ m is a positive integer (mode number), 
and ( _p( E / M ) denotes the principle argument (phase angle) of (n' — n f cos 9) / (n' + n £ cos 9). 

Equation (129)) is the lasing threshold condition for oblique TE and TM modes. Setting t = 9 = 0, it reproduces 
the result given in [§]. Figure [2] shows the plots of g^ E ! M ^> as a function of 9 for a semiconductor gain medium with 
L = 300 /iin and g = 3.4 that is obtained by setting y( E / M ) = 1500nm and using (1251) and (1221) . The following are 
some of the properties that are demonstrated in Fig. [2] 
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FIG. 2: (Color online) Plots of the threshold gain coefficient g as a function of 6 for the TE waves (thick blue curve) and TM 
waves (thin red curve) of wavelength 1500 nm in a slab with L = 300 /im and g = 3.4. The dotted vertical line marks the 
critical angle 9 C 73.6105° at which g takes its maximum value, namely 461.113 cm -1 . 


1. g^ E \ 0) = g( M \ 0), i.e., the slab begins lasing in the TE and TM waves propagating in the normal direction to 
the slab at the same value of the gain. 

2. For 0 < 8 < 90°, g^ E \d) < g( M \d). Therefore, it is more difficult to initiate lasing in the oblique TM modes 
than the oblique TE modes, i.e., it requires more pumping power. 

3. As 9 —> 90°, both g^ E \d) and g^ E \0) tend to zero. This is to be expected, because for 6 —> 90° the optical path 
inside the gain re gion becomes infinitely long. The situation resembles that of spectrally singular whispering 
gallery modes [lj,[l 3 - It corresponds to the infinite radius limit of the latter. 

4. g is a smooth monotonically decreasing function of 8, whereas g <M ' > (8) increases for 8 < 8 C and deceases for 
8 > 8 Cl for some critical angle 8 C = 73.6015°. 

In order to understand the nature of 8 C and the behavior depicted in Fig. (2j we consider typical optically active 
material which satisfies |k| -C g — 1. This allows for a perturbative treatment of (1291) and (1501) . Expanding the 
right-hand side of these equations in a power series in k and ignoring quadratic and higher order terms, we find 


and 


9 {E \0) 

A (£) 


4\/ g 2 — sin 2 9 1 Wg 2 — sin 2 8 + cos0| 
Lg vV - 1 




2k cos 8 
nm{g 2 — 1) 


g {M \0) 

A (M) 


2\J g 1 — sin 2 
Lg 



\Jg 2 — sin 2 8 + g 2 cos 8 
\/g 2 — sin 2 9 — g 2 cos 9 

sin 2 9 4 k cos 9[g 2 — 1 + 008 ( 20 )] 

g 2 Trm(g 2 — 1 )[? 7 2 — 1 + ( g 2 + 1) cos(20)] ’ 



(31) 

(32) 


(33) 

(34) 


where g’ := \J g 2 — sin 2 9 and “ss” stands for the fact that we have ignored quadratic and higher order terms in k. 
Equations (IQTl) and (|33l) provide practically more appealing expressions for the lasing threshold (and CPA) condition 

(Ell). 
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FIG. 3: (Color online) Graphs of the wavelengths as a function of 9 for L = 300 pm, tj = 3.4, re = —10 4 , and 

m = 650, 665, 680. The difference between A (£ ^ and A ■ iM 1 are too small to be visible. 


A closer examination of (1531) shows that its right-hand side diverges whenever sin0 
tan 9 = rj. This suggests that 


v/V 7 ! 2 + lj equivalently 


9 C ~ tan l (rj) =: 6 b , (35) 

where 9 b denotes the Brewster’s angle [29j of the slab in the absence of gain. Indeed, the Brewster’s angle 9 b for the 
sample used in Fig. [2] is 73.6015°. Therefore, our approximate expression (1351) for 9 C is in perfect agreement with the 
exact numerical result. 

Equation (l33l) seems to imply that diverges for 9 = 9 b . We should however recall that this equation is reliable 
only if |re| -C rj — 1. The prediction that takes extremely large values in the vicinity of 9 b conflicts with this 
condition, because g is proportional to re. In fact, it is easy to see from ([29| that g^ is a bounded function for re ^ 0; it 
attains a large but finite maximum value at 9 = 9 C . For the sample used in Figured we find g ( 9 C ) = 461.113 cm' 1 , 
which is about an order of magnitude larger than the attainable gain coefficients for typical high-gain material [30| . 
Therefore, lasing (and CPA) are virtually unattainable for oblique TM modes with 9 in the close vicinity of the 
Brewster’s angle. This is actually quite easy to understand because for this modes the boundaries of the slab have 
very small reflectance (originating from the imaginary part of the complex refractive index.) This in turn reduces the 
internal reflections and the optical path inside the slab, and results in a very large threshold gain. 

In typical situations, i > A. Therefore, according to (El and (l34l) . the mode number m is much greater than 
1. This in turn implies that, for the same mode number, the difference between A^ and A^ M ) is negligibly small. 
Figure [3] shows the plots of as a function of 9 for the same gain medium as in Fig. d and different mode 

numbers m. 


V. ACCOUNTING FOR DISPERSION 

The characterization of spectral singularities that we provide in the preceding section takes the refractive index n 
and the wavenumber k as independent parameters. In this section, we provide a description of the spectral singularities 
for a more realistic situation where we take into account the effect of dispersion. 

Suppose that our slab is made out of a doped host medium of refraction index n o that we can model by a two-level 
atomic system with lower and upper level population densities Ni and N u , resonance frequency wo, damping coefficient 
7 , and the dispersion relation 


2 2 
n — ni — 


Cj 2 — 1 + fy Cj ’ 


( 36 ) 


where Cj := w/wo, 7 := 7 /wo, Cj p := (Ni — N u )e 2 /(m e £oWo), e is electron’s charge, and m e is its mass. We can express 
Cj 2 in terms of the imaginary part reo of n at the resonance wavelength Ao := 2t:c/loq according to w 2 = 2 nnyren + 0 (re 2 ) , 
[a]. Inserting this relation in ( 1361 ) . using ( 1271 ) . and neglecting quadratic and higher order terms in re 0 , we obtain 10] 


re ps 


V « n 0 + K 0 fi(uj) 


«o/2(w) 


(37) 
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FIG. 4: (Color online) Spectral singularities in the TE Modes of a planar slab of thickness L = 300 /im made out of the gain 
medium (1381) for different incidence angles 9. The horizontal line, go = 40 cm -1 , signifies the experimental upper bound on the 
attainable gain. 


where 


/l(w) 


7(1 ~~ k-' 2 ) 

(1 — w 2 ) 2 + 7 2 o ) 2 ’ 


/2(w) 


7 2 6j 

(1 — U ! 2 ) 2 + j 2 6 j 2 ' 


In view of (1261) . we also have kq = — Ao< 7 o/ 47 r. Substituting this equation in (1371) and using the resulting relations 
together with ESI) in (E5l) we can determine the A and go values for the spectral singularities. 

Figures H 0 and [S] show the location of the spectral singularities in the X-g 0 plane for a semiconducting slab of 
thickness L = 300 /xm with the following specifications [30]. 


no = 3.4, Aq = 1500nm, 7 = 0.02, ffo < 40cm 1 . (38) 


Here the last inequality gives the upper bound on the attainable gain for this material. 


VI. SPECTRALLY SINGULAR TE AND TM WAVES 


If we adjust the parameters of our system so that Eq. (1241) holds for a particular wavenumber k , the reflection and 
transmission amplitude diverge. This corresponds to the emergence of purely outgoing waves, 


ao = b 2 = 0 . 

Equations (l24l) and (|39l) together with the boundary conditions listed in Table ITT1 imply 

n — o~ ik z L h 

a 2 — e Oq. 


Furthermore, we can use (El to establish 


e ikz 



(39) 


(40) 


( 41 ) 


Substituting (l39l) (TUT) in the formulas given in Table Q] we can obtain the explicit form of the spectrally singular 
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FIG. 5: (Color online) Spectral singularities in the TM Modes of a planar slab of thickness L = 300 fim made out of the gain 
medium (1381) for different incidence angles 8. The horizontal line, go = 40 cm -1 , signifies the experimental upper bound on the 
attainable gain. 
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Spectrally Singular TE-Fields 

Spectrally Singular TM-Fields 

E x = E z = H y = 0 

E y = b 0 e kxX F + ( n, z) 

tt fro COS 0 ik x x rp (~ \ 

H x = e x E_ (n, z) 

^0 

fro sin 0 ik x x Tp (~ \ 

H z = 7 e F+{n,z) 

6 0 

E y = H x = H z = 0 

E x = boZo cos 9e tkxX F-(-^, z) 

E z = —boZosm.8 e kxX G+(^-,z) 
H v =b 0 e ikxX F+(±,z) 


TABLE III: Components of the spectrally singular TE and TM fields in cartesian coordinates. Here E±(u, z) and G+(u, z) are 
defined by (1121) (1111) . and (1161) holds. 


TE and TM waves. Table Iffl] gives the result of this calculation in terms of the functions, 



e-ik 

zZ 

for 

z < 0, 


F+(u,z) := < 

U+{ u, 

z)/u for 

0 < z < L, 

(42) 


e ik z (z 

-L) 

for 

z> L, 



e ~ikz 

Z 

for 

z < 0, 


li 

U-(VL, 

z) 

for 0 

<z<L , 

(43) 


e ik z {z- 

< 

-L) 

for 

z > L, 



\ e~ if 

Z z z 

for 

z < 0, 


G+(u,z) := 

U+( u, 

z)l 

n for 

>4 

VI 

VI 

o 

(44) 


e ik z (z 

-L) 

for 

z > L, 



U ± ( u, z) := \ [(u - lY- z/L (u + l) z/L ± (u - 1 ) z/L (u + 1 ) 1_Z / L 


which determine the z-dependence of the fields. Note that the k x and k z that appear in Tabic Hill and 
bound to take values that are consistent with Eq. (1241) . In particular, they satisfy 


gifcxX _ 


u+1 yan 9 x /nL , fc ^ / U +l 

,u-iy ’ 6 ^ “ Vu-1 

The transformation properties of F±, G+, and U± under the parity transformation, 

V r 

Z —> L — z. 


z/nL 


(45) 


are 


(46) 


(47) 


provide a direct verification of the P-invariance of the phenomenon of spectral singularities. This is important, because 
it is the basic symmetry responsible for the coherent emission of waves from the left and right boundaries of the slab 
at a spectral singularity ]2(| . We provide a more explicit demonstration of this behavior in the following section where 
we examine the Poynting vector for the spectrally singular waves. 


VII. POYNTING VECTOR AND ENERGY DENSITY AT A SPECTRAL SINGULARITY 

If the slab we consider has the gain coefficient g required for generating a spectral singularity in an oblique TE 
or TM mode, it begins lasing in this mode. In this section we examine the behavior of the time-averaged Poynting 
vector and the energy density inside and outside the slab for the values of the physical parameters of the slab that 
realize such a spectral singularity. The time-averaged Poynting vector and the energy density are respectively given 
by @ 


{S) = \ Re ( ExH *), 

<«> “ (E ■ D* + B ■ H*) = i (e 0 Re[ 3 (^)]| ^| 2 + ^\H\ 2 ) 


(48) 

(49) 
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We calculate these quantities for the spectrally singular modes of our system by substituting the formulas given in 
Table Hill in Eqs. (135)1 and (1351) . This yields 


sin 9 e x — cos 9 e z for z < 0, 

(S(E/M)) = |<^/ M) )|x<( ( 2 ) sin 0 g x + Z (s/M) (z) cos 9 e z for 0 < z < L, 

sin 9 e x + cos 9 e z for z>L , 


<„<*/">> = („(*/")) x / uiE ' M) (') for 0 s 2 s L ’ 

( 1 for 0 < 0 and z > L, 

where the superscripts (E) and (M) refer to the TE and TM waves, and 


1(^)1 := 

X (E) (z) := 


\bo\ 2 


2Z 0 ’ 

U+{n,z) 


\(&\ ■■= 
X^ M \z) := 


Z2\i>oY 


U+(h/n 2 ,z) 


Z (E \z) :=Re 


U+(n, z)U-(n, z)* 


Z^ M \z) := Re 


n 

n 2 U + ( n/n 2 , z)f7(n/n 2 , z) 


Re(n 2 ), 


2 


(50) 


(51) 


/„(E)\ e ol fc ol 2 lJ M \ — ^°N 2 

X^Ollt/ ■ 2 ’ x^out / * 9 

U (E \z) := -4^ | [Re(n 2 ) + sin 2 9] |E/+(tt, z)\ 2 + cos 2 9 |n/7_(n, z)| 2 }, 

U( M \z) := 1 | [|n| 4 + sin 2 #Re(n 2 )] 1 (n/n 2 , z)\" + cos 2 0Re(n 2 ) | n C/_(n/n 2 , z)| 2 j . 

Figure [7] shows plots of ( S ) for a slab of thickness L = 400^m and real refractive index g = 3.4 that admits spectrally 
singular TE and TM modes with 9 = 20° and 


\ {T/E) = 1500.111 nm, 


g (E) = 28.385 cm" 1 , 


g (M) = 32.048 cm" 1 . 


(52) 


Figure [ 8 ] shows the plots of (u)/(u 0 u t) for the spectrally singular modes of Fig. [7] and spectrally singular TE and TM 
modes with 9 = 80° and the following specifications 


\( E ) = 1500.519 nm, g {E) = 5.118 cm" 1 , 
A (m) = 1500.506 nm, g {M) = 68.982 cm" 1 . 


(53) 


These respectively provide the general behavior of ( u) for 9 < 9 C and 9 > 9 C . In particular, for the TM modes with 
9 > 9 C the energy density inside the slab is smaller than that of outside the slab. This is in sharp contrast to TE 
modes and TM modes with 9 < 9 C . The same applies to the magnitude of time-averaged Poynting vector; for TM 
modes with 9 > 9 Cl | (S') | takes much smaller values inside the slab and increases sharply at the boundaries. 

Figures [7] and [5] provide a clear demonstration of the T’-syinmetry of the phenomenon of spectral singularities. It 
is indeed easy to establish the existence of this symmetry using Eqs. (1501) and (1511) . Under the parity transformation 
(1371) . U±{u,z), X^ E / M \ z( E / M \ and transform according to 

U± (u, z)->±U± (u, z ), X^W -> X^W , 

Z (S/M) _ Z (E/M), U (E/M) _^ U (E/M )_ 


These in turn imply that and ( m( b / m )) are P-invariant quantities. 

The following are some of the other notable features of and (u^ E ' M ^) that can be established using the 

analytic expressions given in (l50l) and (l5ll) or seen from Figs. [7] and [8] 

1. (S^) is a continuous function of z. In particular, if we denote the angle between (S^) and the positive z-axis 
by 0 ( b H -) 5 we find that 0 ( E )(O) = tt - 9 and @( e \L) = 9. 
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FIG. 7: (Color online) Graphs of the time-averaged Poynting vectors {S < ' E ' > ) (on the left in blue) and (S^ M ^) (on the right in 
red) for spectral singularities in oblique TE and TM modes with 0 = 20° and specifications (1521) . Here L = 400 /im and r/ = 3.4. 
The gray region represents the slab. 



z (fim) 


0 = 80° 



FIG. 8: (Color online) Graphs of the normalized energy density ( u)/{u 0 ut ) for the spectrally singular TE (solid blue curves) 
and TM (dashed red curves) modes. The left-hand figure corresponds to the configuration depicted in Fig. |T]that has 6 = 20°. 
The right-hand figure corresponds to spectral singularities (53]) with the same values of L and tj, but for 6 = 80°. Again the 
gray area corresponds to the gain region. 


2 . (S' M ') is a continuous function of z except at the boundaries of the slab where its ^-component undergoes 
jumps. This follows from (1501) and the fact that 

0 ) = X ( M )(L) = ^p-^l, 

Z (M) ( 0) = -Z {m \L) = ~ 1. 















13 


In particular, the angle Q^ M \z) between (S < ' M ' > ) and the positive 2 -axis satisfy 


{ 7r 8 for 2 < 0, 

n-e for 2 = 0, 

8 for z = L, 

d for 2 > L, 


6 := tan 1 


tan 8 Re(n 2 ) 

l n l 4 


For the spectrally singular 

3. At 2 = L/2 both ( S {E) ) 
identically; 


TM mode depicted in Fig. [7] where 8 = 20°, we have 8 = 1.8°. 
and ( S( M ' ) ) point along the positive x-axis, except for 8 = 


<S< B >> 


z=L/2 


z—L/2 


= K^out)! |n-n *| 2 sin0e x , 


= \(SZ } )\ 


/Re(n 2 )\ 

ft n 2 

l n 4 ) 

n 2 ft 


sin a e x . 


0 where they vanish 


4. Inside the slab |( I S( B / M ))| and ( u ( e /m )) a ^ a j n their minimum value at 2 = L/2. The latter grows monotonically 
as one approaches the boundary. 

The graphical description of (g( B / M )) and for spectral singularities with different values of physical 

parameters turns out to have general validity. In particular, there is no qualitative difference between oblique spectrally 
singular TE modes with different 8. For the oblique spectrally singular TM modes there are basic differences of 
behavior between the modes with 8 < 8b and those with 8 > 


VIII. CONCLUDING REMARKS 

In this article we have used the recent developments regarding the optical realizations of spectral singularities to 
provide a mathematically precise derivation of the laser threshold condition for the oblique TE and TM modes of an 
infinite planar slab. The standard textbook treatments of this problem rely on the principle of conservation of energy 
(3Cj. Our approach has the advantage of not relying on any physical principles but only on the Maxwell’s equations 
and the definition of a spectral singularity. 

For the TE waves, the essential features of spectral singularities, which determine the lasing threshold and the CPA 
conditions, are not sensitive to the direction of the propagation of the wave. This turns out not to be the case for 
the TM waves. If the incidence (or rather emission) angle 8 is smaller than the Brewster’s angle 8 the spectrally 
singular TM waves behave exactly like spectrally singular TE waves, except that their threshold gain is an increasing 
function of 8. For the TM waves with 8 > 8b the situation is completely different. 

A basic feature of spectral singularities is that they are 'P-invariant. This symmetry is responsible for the fact that 
at a spectral singularity the emitted waves from the left and right boundaries of the slab have identical amplitude 
and phase, i.e., the slab emits coherent waves. We have provided a detailed examination of this phenomenon by 
offering explicit formulas for the energy density and Poynting vector of the oblique spectrally singular TE and TM 
waves. This reveals the curious fact that the energy density of a spectrally singular TM wave with 8 > 8b is smaller 
inside the gain region than the surrounding vacuum. This is the opposite of what happens for spectrally singular 
TM waves with 8 < 8b and the spectrally singular TE waves with arbitrary 8. This observation has counterintuitive 
implications. For example it means that if we use a lossy material and arrange that its attenuation coefficient has 
the same value as the threshold gain coefficient supporting an oblique spectrally singular TM mode with 8 > 8b, the 
slab will function as a CPA having a lower energy density inside, i.e., the waves are absorbed more predominantly 
by its boundary. 
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